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Non-integrability of some Hamiltonians with 
rational potentials 

Primitivo Acosta-Humanez & David Blazquez-Sanz 
7th February 2008 

Abstract 

In this work we compute the families of classical Hamiltonians in two 
degrees of freedom in which the Normal Variational Equation around an 
invariant plane falls in Schrodinger type with polynomial or trigonometri- 
cal potential. We analyze the integrability of Normal Variational Equation 
in Liouvillian sense using the Kovacic's algorithm. We compute all Ga- 
lois groups of Schrodinger type equations with polynomial potential. We 
also introduce a method of algebrization that transforms equations with 
transcendental coefficients in equations with rational coefficients without 
changing essentially the Galoisian structure of the equation. We obtain 
Galoisian obstructions to existence of a rational first integral of the orig- 
inal Hamiltonian via Morales- Ramis theory. 

Key Words: Picard-Vessiot Theory, Hamiltonian systems, Integrability, 

Kovacic's Algorithm. 
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1 Introduction 

In [Mo-Sij . C. Simo and J. Morales- Ruiz find the complete list of two degrees 
of freedom classical hamiltonians, with an invariant plane T — {x 2 =2/2 = 0} 
with Normal Variational Equation of Lame type along generic curves in T. 
In those computations they use systematically differential equations satisfied 
by coefficients of Lame equation. It allow us to develop a method to find the 
families of hamiltonians with invariant plane T, with the property of having NVE 
satisfying certain conditions. It allow us to find the families of hamiltonians 
with generic NVE of type Mathieu, Shrodinger with polynomial potential of 
odd degree, quantum harmonic oscillator, and some other simpler examples. 

Then, we analyze those families of linear differential equations in the con- 
text of Picard-Vessiot theory. For equations with rational coefficients, we use 
Kovacic's Algorithm. In particular we give a complete descripion of Galois 
groups of Schrodinger type equations with polynomial potential (theorem [fjj) . 
For equations with transcendental coefficients, we develop a method of alge- 
brization. We characterize equations that can be algebrized through change of 
variables of certain type (hamiltonian). We prove the following result. 

Algebrization algorithm. The differential equation y — r(t)y is algebriz- 
able through a hamiltonian change of variable x = x(t) if and only if there exists 
f,a such that ^, £ G C(af), where f(x(t)) = r(t), a(x) = 2(H - V{x)) = 
x 2 . Furthermore, the algebraic form of the equation y = r(t)y is 

.. la', f 

y" + --y'--y = o. 

2 a a 



Once we get the complete analysis of linearized equations, we apply a theo- 
rem of Morales-Ramis, obtaining the following results on the non integrability 
of those hamiltonians for generic values of the parameters. 

Non-integrability Results. Hamiltonians, 

2,2 

1. "^-y^ + {X2+ 2x 3Xl yi + x o - AiieS; - X 2 xxxl - \ 3 xlx% + 0(x 1 ,X2)x%, with 
A3 ¥= 0; 
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2. Vl 3 2 + Ao + Q(xi)a;2 + £2)^21 where is a non-constant poly- 
nomial; 

2 + 2 aj 2 a; 2 

3. yi 2 V2 + /l*o + Mi^i + "ir 1 — ^0^2 — X1X1X2 + /3(xi,X2)x2, with w ^ 0. 

4 2/i+2/| 1 1 Ml 1 \i^ 2 x t , ^ 2 z 2 , \ 2 \ „ „2 \ _2_2 , 

4. h A*0 + (A!+2A 2Kl )2 H 8Al 1 8 h A 0»2 ~ 'MZl^ - A2X^X2 + 

(3{x\,X2)x2, with w 7^ 0, and A2 7^ 0. 

where (3(xi,X2) is any analytical function around T = {X2 = y% = 0}, do no 
admit any additional rational first integral. 

Corollary. Every integrable (by rational functions) polynomial potentials 
with invariant plane T = {X2 = g/2 = 0} can be written in the following form 

V — Qi(xi, x 2 )x\ + Xix\ + Ao, Ao,AigC. 

Note that for (3(x\,X2) polynomial and some values of the parameters we 
fall in the case of homogeneous polynomial potentials. Those cases had been 
analyzed deeply in |Ma-Prl INa-Yoj . 



1.1 Picard-Vessiot theory 

The Picard-Vessiot theory is the Galois theory of linear differential equations. 
In the classical Galois theory, the main object is a group of permutations of 
the roots, while in the Picard-Vessiot theory is a linear algebraic group. For 
polynomial equations we want a solution in terms of radicals. From classical 
Galois theory it is well known that this is possible if and only if the Galois group 
is solvable. 

An analogous situation holds for linear homogeneous differential equations. 
For more details see |Va-Si| . The following definition is true for matrices n x n, 
but for simplicity we are restricting to matrices 2x2. 

Definition 1 An algebraic group of matrices 2x2 is a subgroup G C GL(2, C), 
defined by algebraic equations in its matrix elements. That is, there exists a set 
of polynomials 

{P l {xn,Xi2 1 X21 1 X 2 2)}ieI, 

such that 

11 X12 gG ^ Vie I,Pi{xn,Xl2,X21,X22) = 0- 

X21 X22 ) 

In this case we say that G is an algebraic manifold provided of an structure of 
group. In the remainder of this paper we only work, as particular case, with 
linear differential equations of second order 

y" + ay' + by — 0, a, b G C(x). 
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Suppose that yi , y 2 is a fundamental system of solutions of the differen- 
tial equation. This means that y\ , yi are linearly independent over C and 
every solution is a linear combination of these two. Let L = C(x){yi,y2) — 
C(x)(yi, 7/2, Vxi that is the smallest differential field containing to C(x) and 
{yi,V2}- 

Definition 2 (Differential Galois Group) The group of all differential au- 
tomorphisms of L over C(x) is called the Galois group of L over C(x) and 
denoted by Gal(L/C(x)) or also by GclIqm- This means that for a: L — ► L, 
o~(a!) = er'(a) and Va € C(x), a (a) = a. 

If a £ Gal(L/C(x)) then ctj/i, crj/2 is another fundamental system of solutions 
of the linear differential equation. Hence there exists a matrix 

A =(c h d)^ GL ^ 

such that 

a ( Vl \ = ( ayx \ = A ( y A . 
\V2j \o-y2j \y2j 

This defines a faithful representation Gal(L/C(x)) — > GL(2, C) and it is possible 
to consider Gal(L/C(x)) as a subgroup of GL(2, C). It depends on the choice 
of fundamental system yi, y 2 , but only up to conjugate. 

One of the fundamental results of the Picard-Vessiot theory is the following 
theorem. 



Theorem 1 The Galois group G = Gal(L/C(x)) is an algebraic subgroup of 
Gi(2,C). 

Now we are interested in the differential equation 

e = r£, reC(x). (1) 

We recall that equation can be obtained from the general second order 
linear differential equation 

y" + ay' + by = 0, a,beC(x), 

through the change of variable 

y = e 2 J a £ r = 1 b. 

4 2 

On the other hand, through the change of variable v — we get the 
associated Rlccatti equation to equation (|T|) 

v' = r — v 2 , v = (2) 
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For the differential equation Q, G = Gal(G/C(x)) is an algebraic subgroup 
SL(2,C). 

Recall that an algebraic group G has a unique connected normal algebraic 
subgroup G° of finite index. This means that the identity component G° is the 
biggest connected algebraic subgroup of G containing the identity. 

Definition 3 Let F be a differential extension ofC(x), and let be rj solution of 
the differential equation 

y" + ay' + by = 0, a,b € F 

1. X] is algebraic over F if r\ satisfies a polynomial equation with coefficients 
in F , i.e. rj is an algebraic function of one variable. 

2. rj is primitive over F if rj' £ F , i.e. f] = J f for some f € F. 

3. r) is exponential over F if ' rf /rj G F , i.e. rj — e-> f for some f £ F. 

Definition 4 A solution rj of the previous differential equation is said to be 
Liouvillian over F if there is a tower of differential fields 

F = F cF 1 c...cF m = L, 

with n £ L and for each i = l,...,m, Fi — Fi-\(jji) with rji either algebraic, 
primitive, or exponential over Fi—\. In this case we say that the differential 
equation is integrable. 

Thus a Liouvillian solution is built up using algebraic functions, integrals and 
exponentials. In the case F = C(x) we get, for instance logarithmic, trigono- 
metric functions, but not special functions such that the Airy functions. 

We recall that a group G is called solvable if and only if there exists a chain 
of normal subgroups 

e = G < Gi < . . . < G n = G 
such that the quotient GijGj is abelian for all n > i > j > 0. 

Theorem 2 The equation Q is integrable (has Liouvillian solutions) if and 
only if for G = Gal(L/C(x)), the identity component G° is solvable. 

Using the theorem [3 Kovacic in 1986 introduced an algorithm to solve the 
differential (|TJ and show that JTJ is integrable if and only if the solution of the 
equation J2J) is a rational function (case 1), is a root of polynomial of degree two 
(case 2) or is a root of polynomial of degree 4, 6, or 12 (case 3) (see Appendix 
A). Based in the Kovacic's algorithm we have the following key result. 

The Galois group of differential equation with r = Qk{x) a polynomial of 
degree k > is a non-abelian connected group. 

For complete result, details and proof see Appendix A. 



5 



1.2 Morales-Ramis theory 



Morales- Ramis theory jMo-Rall IMo-Ra2"] . see also |Mol| . relates the integrabil- 
ity of hamiltonian systems to the integrability of linear differential equations. In 
this approach analyze the linearization (variational equations) of hamiltonian 
systems along some known particular solution. If the hamiltonian system is in- 
tegrable, then we expect that the linearized equation has good properties in the 
sense of Picard-Vessiot theory. Exactly, for integrable hamiltonian systems, the 
Galois group of the linearized equation must be virtually abelian. It gives us the 
best non-integrability criterion known for hamiltonian systems. This approach 
has been extended to higher order variational equations in [M*-R-S . 



1.2.1 Integrability of hamiltonian systems 

A symplectic manifold (real or complex), M 2n is a 2n-dimensional manifold, 
provided with a closed 2-form u. This closed 2-form gives us a natural isomor- 
phism between vector bundles, b: TM — > T*M. Given a function H on M, 
there is an unique vector field Xh such that, 

b(X H ) = dH 

this is the hamiltonian vector field of H. Furthermore, it gives an structure of 
Poisson algebra in the ring of differentiable functions of M% n by defining: 

{H, F} = X H F. 

We say that H and F are in involution if {H, F} = 0. From our definition, it is 
obvious that F is a first integral of Xh if and only if H and F are in involution. 
In particular H is always a first integral of Xh- Moreover, if H and F are in 
involution, then their flows commute. 

The equations of the flow of Xh, in a system of canonical coordinates, 
pi, . . . ,p n ,qi, ■ ■ ■ ,q n (i-e. such that uj 2 = Yh=i P l A are written 

OH , ^ x OH , 

<i=-^(={H,q}), p = -—(={H,p}), 

and they are known as Hamilton equations. 



Theorem 3 (Liouville- Arnold) Let Xh be a hamiltonian defined on a real 
symplectic manifold Mi n . Assume that there are n functionally independent first 
integrals Fi,...,F n in involution. Let M a be a non-singular ( i. e. dF\ , . . . , dF n 
are independent on very point of M a ) level manifold, 

M a = {p: Fi{p) = ai, ... , F n {p) = a n }. 

Then, 

1. If M a is compact and connected, then it is a torus M a ~ R"/Z™. 
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2. In a neighborhood of the torus M a there are functions I\, 4>\, . . . , 4> n 

such that 

n 

L02 = dh A dfa, 

and {H, Ij} = for j = 1, . . . , n. 

From now on, we will consider C 2n as a complex symplectic manifold. 
Lioville- Arnold theorem gives us a notion of integrability for Hamiltonian sys- 
tems. A hamiltonian H in C 2 " is called integrable in the sense of Liouville if 
there exist n independent first integrals of Xh in involution. We will say that 
H in integrable by terms of rational functions if we can find those first integrals 
between rational functions, or whatever. 

1.2.2 Variational equations 

We want to relate integrability of hamiltonian systems with Picard-Vessiot the- 
ory. We deal with non-linear hamiltonian systems. But, given a hamiltonian 
H in C 2 ™, and T an integral curve of Xh, we can consider the first variational 
equation (VE), as 

£xA = o, 

in which the linear equation induced in the tangent bundle (£ represents a vector 
field supported on T). 

Let r be parameterized by 7: 1 1— > {x(t), y(t)) in such way that 

dxi dH dy.i dH 
dt dyi ' dt dxi 

Then the VE along T is the linear system, 

ii\ = ( sSftw*)) s£ft(7(*)) 

ft'/ \-Bg§7(7(t)) -*gfe(7(*)) 
From the definition of Lie derivative, comes that 

is a solution of the VE. We can use it to reduce, using some generalization of 
DAlambert method |Mo-Rall IMo-Ra2] (see also |Molj ) our equation, obtaining 
the Normal Variational Equation (NVE), which is a linear system of rank 2(n — 
1). In the case of 2-degrees of freedom hamiltonian systems, those NVE can be 
seen as second order linear homogeneous differential equation. 

1.2.3 Non-integrability tools 

Morales- Ramis theory, concerns several results that relate the existence of first 
integrals of H with the Galois group of the variational equations. Along this 
paper we will use systematically the following one: 
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Theorem 4 ([Mo-Ral ) Let H be a Hamiltonian in C 2n , and 7 a particular 
solution such that the NVE has irregular singularities at the points of 7 at 
infinity. Then, if H is completely integrable by terms of rational functions, then 
the Identity component of Galois Group of the NVE is abelian. 

Remark 1 Here, the field of coefficients of the NVE is the field meromorphic 
functions on 7. 

2 Some results on linear differential equations 

2.1 Algebrization of Linear Differential Equations 

For some differential equations it is useful, if is possible, to replace the original 
differential equation over a compact Riemann surface, by a new differential 
equation over the Riemann sphere P 1 (i.e., with rational coefficients) by a change 
of the independent variable. This equation on P 1 is called the algebraic form 
or algebrization of the original equation. Kovacic's algorithm can be applied 
over the algebraic form to solve the original equation. In a more general way we 
will consider the effect of a finite ramified covering on the Galois group of the 
original differential equation. In Moll lMo-Ralj the following theorem is given. 

Theorem 5 (Morales- Ramis |MolL iMo-Ralj ) Let X be a (connected) Rie- 
mann surface. Let f : X' — > X) be a finite ramified covering of X by a Riemann 
surface X' . Let V be a meromorphic connection on X . We set V' = /*V. Then 
we have a natural infective homomorphism Gal(S7') — ► Gal(V) of differential 
Galois groups which induces an isomorphism between their Lie algebras. 

Remark 2 We observe that, in terms of the differential Galois groups, this 
theorem means that the identity component of the differential Galois group is 
invariant by the covering. In other words, if the original differential equation 
over the Riemann surface X' is transformed by a change of the independent vari- 
able in a differential equation over the Riemann surface X , then both equations 
have the same identity component of the differential Galois group. 

Proposition 1 (Change of the independent variable) Let us consider the 
following equation, with coefficients in C(x): 

y" + a(x)y' + b(x)y = 0, y' = ^ (3) 

and C(x) <^-> L the corresponding Picard-Vessiot extension. Let (K,S) be a 
differential field with C as field of constants. Let £ £ K be a non- constat. 
Then, by the change of variable x — £, © is transformed in, 

y+ ] y + m(i) 2 y = o, * = ?>z. (4) 
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Let, Kq C K be the smallest differential field containing £ and C. Then Q *s a 
differential equation with coefficients in Kq. Let Kq Lq be the corresponding 
Picard-Vessiot extension. Assume that 

C(x) ^K , x^Z 

is an algebraic extension, then 

Gal{L a /K ) Q = Gal(L/C(x))°. 

Proof. By the chain rule we have 

— = ls, 

dx £ ' 



and, 



dx* (£)2 (£)3 



now, changing y', y" in J2J and making monic this equation we have Q 

In the same way we can obtain Q through Q . 

If is an algebraic extension of C(x), then we can identify Kq with the 
ring of meromorphic functions on a compact X Riemann surface, and £, is a 
finite ramified covering of the Riemann sphere, 

X^S 1 . 

If we consider V the meromorphic connection in S 1 induced by equation 
then £*(V), is the meromorphic connection induced by equation J3J in X, and 
finally by theorem we conclude. □ 

Recently Manuel Bronstein in |B] has implemented an algorithm to solve 
differential equation over C(t, e' ^ t ') without algebrize the equation. As imme- 
diate consequence of the proposition J2J we have the following corollary. 

Corollary 1 (Linear differential equation over C(t,eS Let be f € C(t), 

the differential equation 



y-(f(t) + i| - f(t)eJ /Wo (ef '<«>) ) (V «*)))" 6 /W) „ 



by the change x = e-> /W is algebrizable and its algebraic form is given by 

y" + a(x)y' + b(x)y = 0. 



-0, 
(5) 
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Remark 3 In this corollary, we have the following cases. 

1. / — rip for a rational function h, n € Z + , we have the trivial case, both 
equations are over the Riemann sphere and they have the same differential 
field, so that the equation J^J do not need be algebrized. 

2. / = — h-, for a rational function h, n € Z + , the equation (J5J is defined over 
an algebraic extension of C(t) and so that this equation is not necessarily 
over the Riemann sphere. 

3- / 7^ l\i f° r an Y rational function ft,, g G Q, the equation (JjjJ is defined 
over an transcendental extension of C(t) and so that the this equation is 
not over de Riemann sphere. 

In first and second case, we can apply proposition^] taking Kq = C(i, e-> ^^), 
so that the identity component of the algebrized equation es conserved. The 
conservation of Galois group in the third case, requires further analysis, and will 
not be discussed here. We just remark that the Galois group corresponding to 
original equation is a subgroup of the Galois group of the algebrized equation. 

Definition 5 (Regular and irregular singularity) The point x — xq is called 
regular singular point (or regular singularity) of the equation JSJ) if and only if 
x — xq is not an ordinary point and 

(x-xo)a(x), (x — xo) 2 b(x), 

are both analytic in x = Xq. If X = Xq is not a regular singularity, then is an 
irregular singularity. 

Remark 4 (Change to infinity) To study asymptotic behaviors (or the point 
x = oo) in @ we can take x(t) — j and to analyze in the behavior in t = 0. 
That is, to study the behavior of the equation J2J in x = oo we should study 
the behavior in t — of the equation 

2 / 1 \ . 1 . /I 



t ^j a {-t))y + ^ b {-t) y=0 - ^ 

In this way, by the definitional we say that x = oo is a regular singularity of 
the equation if and only if t = is a regular singularity of the equation © . 

To algebrize second order linear differential equations is easier when the term in 
y is absent and the change of variable is hamiltonian, that is, RLDE y = r(t)y. 

Definition 6 (Hamiltonian change of variable) A change of variable x = 
x(t) is called hamiltonian if and only if (x(t),x(t)) is a solution curve of the 
autonomous hamiltonian system 

H = H(x,y) = ^ + V(x). 
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Assume that we algebrize equation (0J through a hamiltonian change of 
variables, x — £(i). Then, K$ = C(£, £,...), but, we have the algebraic relation, 

(0 2 = 2h-2V@, h = H(S,£)eC, 

so that Kq — C(£, £) is an algebraic extension of C(x). We can apply proposition 
2] and then the identity component of the Galois group is conserved. 

Proposition 2 (Algebrization algorithm) The differential equation 

jj = r(t)y 

is algebrizable through a hamiltonian change of variable x — x(t) if and only if 
there exists /, a such that 

a' f 

— , - G C(x), where f(x(t)) = r(t), a(x) = 2(H - V(x)) = x 2 . 
a a 

Furthermore, the algebraic form of the equation y = r(t)y is 

y"+\-y'- f ~y = u. (J) 

2 a a 

Proof. Because x = x(t) is a hamiltonian change of variable for the dif- 
ferential equation y — r(t)y so that x — y, y = x = —V'(x) and there exists 
/, a such that y — f(x(t))y and x 2 = 2(H — V(x)) = a(x). By the proposi- 
tion have —f(x) = b(x)x 2 and a(x)x — xjx — 0, therefore a(x) — x/x 2 
and b(x) = —f(x)/a(x). In this way y = r(t)y is algebrizable if and only if 
a(x),b(x) G C(x). As 2{H - V(x)) = a{x) then a'(x) = -2V'{x) = 2x, and 
therefore a(x) = f ■ In this way, we obtain the equation 0. □ 

As consequence of the proposition [21 we have the following result. 

Corollary 2 Let be r(t) = g(%\, ■ ■ ■ ,x n ), where Xi = e Xit , A, G C. The differ- 
ential equation y = r(t)y is algebrizable if and only if 

^eQ, l<i^j<n, geC(x). 

Furthermore, we have A 2 ; = CiX, where A € C and Cj G Q and one change of 
variable is 

x = e~ , where Ci — — , gcd(pi,qi) = 1 and q = mcm{q\, ■ ■ ■ ,q n )- 

Remark 5 (Using the algebrization algorithm) To algebrize the differen- 
tial equation y = r(t)y it should keep in mind the following steps. 

Step 1 Find a hamiltonian change of variable x = x{t). 
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Step 2 Find / and a such that r(t) = f{x(t)) and {x(t)) 2 = a(x(t)). 
Step 3 Write f(x) and a(x). 

Step 4 Verify if f(x)/a(x) g C(x) and a'(x)/a(x) € C(x) to see if the RLDE 
is algebrizable or not. 

Step 5 If the RLDE is algebrizable, write the algebraic form of the original 
equation such as follows 

,, la', f 

y "+ ~-v '--y = o. 

2 a a 

When we have algebrized second order linear differential equation we study its 
integrability and its Galois groups. 

Remark 6 (Monic Polynomials) Let be the reduced linear differential equa- 
tion (RLDE) 

y=[^2 c kt k )y, cfc e c, fc = i,...,n. 

By the algebrization algorithm we can take x = [it, [i S C, so that 

x = fi, a(x) = i? , a'(x) = and f(x) = z2 Ck \~) ' Cfc ' ^ e C 

fc=o 

Now, by the equation JTJ) the new differential equation is 

«^ = (e (^) **) * c ^ eC - 

In general, for /i — "+^/c^ we can obtain the equation 

V"=[x n +1^!/, <**=(^), fc = 0,.-,n-l. 

Furthermore, we can observe, by definition|Sland equation ©, that the point 
at oo is an irregular singularity for the differential equation with non-constant 
polynomial coefficients because using equation © we can see that zero is not an 
ordinary point and neither is a regular singularity for the differential equation. 

Remark 7 (Extended Mathieu) The Extended Mathieu differential equa- 
tion is 

y = (a + bsini + ccost)y, (8) 

in particular, when b = or c = and \a\ + \b\ ^ 0, we have the so called 
Mathieu equation. Applying the corollary [21 and the steps of the algorithm we 
have x = e lt , x = ix, therefore 

(b + c)x 2 + lax + c - b 2 , 

fix) = — , a{x) = -x\ a'(x) = -2x, 
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so that the algebraic form of the equation (JSJ is 

1 , (b + c)x 2 + 2ax + c - b 

-V +" " 7T 

x Ix 

Making the change x = 1/z in the equation @ we obtain 



y + -y + — ^ y = o, (9) 



Ki)K (c -^ + 2 2 r <t+c) >=°- <"» 

We can observe, by definition [5] and equation ©, that 2; = is an irregular 
singularity for the equation (|10fl and therefore x = 00 is an irregular singularity 
for the equations @ and (JHJ). 

Now, we compute the Galois group and the integrability in equation (0 . So 
that, the RLDE is given by 

£/, = _ ^ (b + c)x 2 + (2a+l)x + c-b ^ ^ (n) 

Applying the Kovacic's algorithm, see Appendix A, we can see that for b =/= —c 
this equation falls in case 2: (03, 003), Eo = {3}, = {1} and so that D = 
because m = — 1 ^ Z + . In this way we have that the equation is not 
integrable, the Galois Group is the connected group SL(2, C), and finally, by 
the theorem the identity component of the Galois group for equation (JSJ) is 
exactly SL(2, C) that is not an abelian group. In the same way, for b = —c we 
have the equations 

y={a-be u )y, y" + -y' + y = 0, 

x lx* 

in which 00 continues being an irregular singularity. Now, Its RLDE is given by 

„ = _ n2a+l)x-2b \ 
\ 2x 3 J 

applying Kovacic's algorithm we can see that this equation falls in case 2: 
(c 3 ,oo 2 ), E = {3}, Eoo = {0,2,4}, so that D = because l/2(eoo - e ) $ Z+. 
This means that the Galois group continues being SL(2, C). Using this result, 
taking e instead of i, we can say that in the case of harmonic oscillator with 
exponential waste 

y = {a + be~ et )y, e > 0, 

the point at 00 is an irregular singularity and the identity component of the 
Galois group is SL{2, C). 

2.2 Galois groups of Schrodinger equations with polyno- 
mial potential 

Let be the reduced linear differential equation (RLDE) 



\k=0 
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By the remark EH through the change of variable x = fit, fi = n+ -^/ck, this 
equation become in 

£" = P n (x)^, P n {x) is a monic polynomial of degree n. 

Kovacic in |Ko| remarked that for n = 2k + 1 there are not liouvillian so- 
lutions and therefore the Galois group of the RLDE is SL(2,C). Using the 
Kovacic's algorithm (see Appendix), we can see that for n = 2k the equa- 
tion falls in case 1, specifically in cq (because has not poles) and 003 (because 
= ~2k), that is {00,003}. 

Lemma 1 (Completing Squares) Every monic polynomial of degree even 
can be written in one only way completing squares, that is 

2n-l / n-1 \ 2 n-1 

Q2n{x) = X 2n + 1kX k = ix U + ^ + ^ • 

k=0 \ k=0 ) k=0 

Proof. Firstly, we can see that 

(n— 1 \ 2 Tl—l n — In — 1 

x n + a kX k = x 2n + 2 a k x n+k + a ^x k +^ , 

fe=0 / fe=0 fc=0 j=0 

so that by indeterminate coefficients we have 

_ <?2n-l _ l2n-2 — Q n -1 _ <?2n-3 ~ 2a„_ia„_2 

f*n— 1 ~ ) ^n — 2 ~ 5 ^n— 3 ~ ; ) 



q n — 2a\a n -\ — 2a2a„_2 — ■ ■ ■ , 2 u o 

a = , b = q -a , b 1 = q 1 - 2a a 1: 

b n -i = 9)1-1 — 2aoa„_i — 2a\a n -i — ■ ■ ■ , 
therefore, we has proven the lemma. □ 

By case 1, {co,oo 3 }, of the kovacic's algorithm, remark ^] [15] (See Ap- 
pendix), remark Eland by lemma have proven the following theorem. 

Theorem 6 (Galois groups in polynomial case) Let us consider the equa- 
tion, 

with Q(x) a polynomial of degree k > 0. Then, it falls in one of the following 
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1. k — In is even, ±6 n _j — n = 2m, m € Z+, and there exist a monic 
polynomial P m of degree m satisfying 

(n-l \ / n-2 n-l \ 

x n + *kx k P' m + nx n - x + Y,( k + 1 H+iX k + Y, b ^ k Pm = 0, 
fc=0 / V k=0 k=0 I 

the solutions are given by 

£l = p m e "+ 1 + ^<==o fc+i , 
6 = P m e "+ 1 + ^*=° 



=0 k+1 

,-,0 2 ( "n+i +£k=0 

and f/ie Galois group is C* X C (non-abelian, resoluble, connected group). 
2. The equation has not liouvillian solutions and its Galois group is SL(2, C). 

Remark 8 (Quadratic case) consider the case where r is a polynomial of 
degree two: 

y = (At 2 + Bt + C) y. 

There are no poles and the order at oo, or^, is —2, so we need to follow case 1 
in {00,003} of the algorithm. Now, by remarkEJand by lemmaG]we have 

y"= {(x + a) 2 + b) y. 

We find that 

[y/r]oo = x + a 
a± =i(±6-l) 
m = o;+ or a~. 

If b is not an odd integer then m cannot be an integer so case 1 cannot hold so 
the RLDE in x has no Liouvillian solutions. If b is an odd integer than we can 
complete steps 2 and 3 and actually (only) find a solution it which is 

y = P m e^+ ax . 
In particular, the quantum harmonic oscillator 

y" = (x 2 + X)y 

is integrable when A is an odd integer and the only one solution obtained by 
means of kovacic's algorithm is given by 

y = H m e 2 , 

where H m denotes the classical Hermite's polynomials. 

For another approach to this problem see Vidunas [Vij and Zoladek in |Zo| . 
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3 Determining families of Hamiltonians with spe- 
cific NVE 

Let us consider a two degrees of freedom classical hamiltonian, 

V is the potential function, and it is assumed to be analytical in some open subset 
of C 2 . The evolution of the system is determined by Hamilton equations: 



X\ = 2/i, x 2 = 2/2, 



dV . dV 
dx\ ' ^ 2 dx2 ' 

Let us assume that the plane T = {xi = 0, 2/2 = 0} is an invariant manifold 
of the hamiltonian. We keep in mind that the family of integral curves lying 
on r is parameterized by the energy h — H\t, but we do not need to use it 
explicitly We are interested in studying the linear approximation of the system 
near T. Since T is an invariant manifold, we have 



dV 

dx 2 



= 0, 



so that the Normal Variational Equation for a particular solution 

1 1 ► j(t) = (xi(t),yi = ±i(t),X2 — 0,2/2 = 0), 

is written, 



£ = v, v 



d 2 v 



Let us define, 



(xi) = V(xi,0), a(xi) 



8 2 V 
dx\ 



(a:i,0), 



and then we write the second order Taylor series in x 2 for V, obtaining the 
following expression for H 



H 



vl + vi 



} +<p(xi)-a(x 1 )-^-+(3(x 1 ,x 2 )xl, 



(12) 



which is the general form of a classical a Hamiltonian, with invariant plane 
T.The NVE associated to any integral curve lying on T is, 



(13) 
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3.1 General Method 

We are interested in compute Hamiltonians of the family l|12f> . such that its 
NVE belongs to a specific family of Linear Differential Equations. Then 
we can apply our results about the integrability of this LDE, and Morales-Ramis 
theorem to obtain information about the non-integrability of such Hamiltonians. 

From now on, we will write a(t) = a(xi(t)), for a generic curve 7 lying on 
r, parameterized by t. Then, the NVE is written 

£ = a(i)e (14) 

Problem. Assume that a(t) is a root of a given differential polynomial 
Q(a, a, a, ...) S C[a, d, a, . . .]. We want to comvute all hamiltonians in l(T^ 
satisfying such a condition. 

In this section we give a method to compute, for any given Q(a, d, . . .), the 
family of classical hamiltonians with invariant plane T such that, for any integral 
curve lying on T, the coefficient a(t) of the NVE satisfies, 

Q(a,d,d, . . .) = 0, (15) 

by solving certain differential equations. This method lies under the calculus 
done by J. Morales and C. Simo in |Mo-Sij . 

We should notice that, for a generic integral curve ^(t) — (xi(t),yi — ii(t)) 
lying on T, equation lfH|) depends only of the values of functions a, and <f). It 
depends of a(xi), since a{t) = a(x\{t)). We observe that the curve j(t) is a 
solution of the restricted Hamiltonian, 

h='^ + 4>{x l ) (16) 
whose associated Hamiltonian vector field is, 

X — y — d 4> d 
1 ' dxi dx\ dyi ' 

thus xi(t) is a solution of the differential equation, 'x\ = and then, the 

relation of x±(t) is given by <f>. 

Since j(t) is an integral curve of Xh, for any function f(x±, y\) defined in T 
we have 

where 7* denote the usual pullback of functions. Then, using a(t) — 7* (a), we 
have for each k > 0, 

~dt h 
so that, 

Q(a, a,a,...) = Q( 7 », 7* (X h a), -f(Xfa), . . .). 

There is an integral curve of the Hamiltonian passing through each point of 
F, so that we have proven the following. 



-JW=r(X k h a), (18) 



17 



Proposition 3 Let H be a Hamiltonian of the family l|12l) . and Q(a, a, d, . . .) a 
differential polynomial with constants coefficients. Then, for each integral curve 
lying on T, the coefficient a(t) of the NVE (I14|) verifies Q(a,a,d, ...,)= 0, if 
and only if the function 

Q(xi,yi) = Q(a,X h a,X\a, . . .), 

vanish on T . 



Remark 9 In fact the NVE of a integral curve depends on the parameteriza- 
tion. Our criterion does not depend on any choice of parameterization of the 
integral curves. This is simple, the NVE corresponding to different paramctcr- 
izations of the same integral curve are related by a translation of time t. We 
need just observe that a polynomial Q(a,a,d, . . .) with constant coefficients is 
invariant of the group by translations of time. Then if the coefficient a(t) of 
the NVE (|14fl for certain parameterization of an integral curve ^(t) satisfied 
{Q = 0}, then it is also satisfied for any other right parameterization of the 
curve. 

Next, we will see that Q(x%, y%) is a polynomial in y\ and its coefficients are 
differential polynomials in a, <f>. If we write down the expressions for successive 
Lie derivatives of a, we obtain 

da , „, 

X h a = Vl — , 19 
dx i 

r,d 2 a dd> da ,„„. 



d 3 a ( d ( dd> da 



xla = y "dxj - m [dk lj£ £r ) + 2 £t m ] (21) 



^ 4 4 d 4 a; 2 f d ( d / d(f> da \ ^ d<f> d 2 a\ ^d 3 a d<fi 
h 1 dx\ 1 \dx\ \dx\ \dx\dx\) dx\dx\) dx\dx\ 



'tf> da ^ 
dx\ \dx\ dx\ ) dx\ dx\ j dx 



In general form we have, 



dXl)a dd> dXH 



A-h ct — yi— - — , {Z6) 

dxi dxi oy\ 

it inductively follows that they all are polynomial in y\ with coefficients differ- 
ential polynomials in a, <j). If we write it down explicitly, 

Xla= (24) 

n>fc>0 
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we can see that the coefficients E„ t k{a,(f>) <E C o^^^tf,^ 
following recurrence law, 

E n+ i. k {a,4>) = -^-E nk -i(a,(j>) - (ft + V)E n ,k+i{a,<t>)^— 
ax i dx\ 



satisfies the 



(25) 



with initial conditions, 



da 



ax i 



(26) 



Remark 10 The recurrence law ()25|l and initial conditions Ij26(l determine the 
coefficients E ni k(ct, <j>). We can compute the value of some of them easily: 



• E n<n {a,4>) = for all n> 1. 



E n ,k(a, 



if n — ft is odd, or ft < 0, or ft > n. 



3.2 Some Examples 

Here we compute families of hamiltonians (|12f> such give rise to specific NVE. 
Although, in order to this computations, we need to solve Polynomial Differ- 
ential Equations, we will see that we can deal with this in a branch of cases. 
Particularly, when Q is a Differential Linear Operator, we will obtain equations 
that involve products of few Linear Differential Operator. 



Example 1 Harmonic oscillator equation is 



(27) 



with Co constant. Then, a hamiltonian of type l|12|) gives such NVE if a = Q. 
Looking at formula (|19|) . it follows that = 0, so that a is a constant. We 
conclude that the general form of a Hamiltonian l|12fl which give rise to NVE of 
the type Q2J is, 

H= Vl * V2 +eP(x 1 ) + X xl+P(x 1 ,x 2 )x 3 2 , 
being Aq a constant, and 0, /? arbitrary analytical functions. 



Example 2 In |Au| . M. Audin notice that the hamiltonian, 



Vx + vl 



gives an example of a simple non-integrable classical hamiltonian, since its NVE 
along any integral curve in L is an Airy equation. Here we compute the family 
of classical hamiltonians that have NVE of type Airy for integral curves lying 
on r. General form of Airy equation is 



£ = (co + cit)£ 



(28) 
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with Co, ci ^ two constants. If follows that a hamiltonian gives rise to NVE 
of this type if a = 0, and d ^ 0. The equation d — gives, by proposition [3] as 
we see in formula (|2UJ), the following system: 



d 2 a dcp da 

dx\ dx\ dx\ 



2-0, ^—=0- (29) 



It split in two independent systems, 




(30) 



Solutions of the first one fall into the previous case of harmonic oscillator. Then, 
taking the general solution of the second system, we conclude that the general 
form of a classical hamiltonian of type (|12() with Airy NVE is: 

H = Vl ^ V ' 2 + A + Xixj + \ 2 x x x\ + /3(xi , x 2 )xl, (31) 

with A2 7^ 0. 

3.2.1 NVE quantum harmonic oscillator 

,3 ,2 

Let us consider now equations with = 0, and %y ^ 0, it is 

|= (c + Cl t + c 2 t 2 )£ (32) 

with C2 =/= 0. Those equation can be reduced to a quantum harmonic oscillator 
equation by an affine change of t, and its integrability has been studied using 
Kovacic's Algorithm. Using proposition and formula (|21|l . we obtain the 
following system of differential equations for a and <fi: 



d 3 a d a d 2 <b d 2 a d<j> 

dx\ dx\ dx\ dx\ dx\ 

General solution of the first equation is 

Al A2 A3 2 

*=- + -x x + -x x , 

and substituting it into the second equation we obtain a Linear Differential 
Equation for </>, 

d 2 </> | 3 2A 3 d<f> =Q 
dx\ A2 + 2A3X1 dxi 

this equation is integrated by two quadratures, and its general solution is 



(A 2 + 2A 3 xi) 
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We conclude that the general formula for hamiltonians of type (|12|) with NVE 
l|32l) for any integral curve lying on T is 



H = vl + vl 



' + Ao — \ix% — A2X1X2 — \ax\x% + (3(xi, X2)x%, (33) 



(A 2 + 2A 3 a;i) 2 
with A3 =/= 0. 



Remark 11 Formula (!.').' jfl is the first example, in this paper, in which we find 
non- linear dynamics in the invariant plane T. Notice that this dynamic is con- 
tinuously deformed to linear dynamics when A4 tends to zero. In general case, 
for a fixed energy h, we have the general integral of the equation: 

8\lh 2 (t - t ) 2 = h{\ 2 + 2A 3 xi) 2 - A 4 . 

3.2.2 NVE with polynomial coefficient a(t) of odd degree 

Let us consider for n > the following differential polynomial, 

d m a 



i(a,a, . . .) 



dt n 



It is obvious that a(t) is polynomial of degree n if and only if Q n (o-, a, . . .) 7^ 
and Q„+i(a,d, . . .) = 0. 

Looking a proposition we see that a hamiltonian p2(l has NVE along a 
generic integral curves lying on T, 

£=Pn(t)Z, (34) 

with P n (t) of degree n, if and only if XJ^a 7^ and X^ +l a vanish on T. Let us 
remember expression ()24ll . X^ +1 a vanish in V if and only if (a, cf)) is a solution 
of the differential system, 

-Rn+l = {-Eri+1,0 («, 4>) = 0) • • • ) E n+ x, n +i{ot, 4>) =0}. 

Using the recurrence law (|25|l defining the differential polynomials E n ,k(a, (f) we 
are going to compute the family of hamiltonians giving rise to a(t) polynomial 
of odd degree. 

Lemma 2 Let (a, <f>) be a solution of R2 m - Then, if 7^ 0, then (a,(f>) is a 

solution fl/]?2m-l- 

PROOF. By remark UHlE^ — 1 oi.(a. 4>) — for all m - 1 > k > 0. Then let 
us proof that -E2m-i,2fc+i(a, 4>) for all m — 2 > k > 0. 
In the first step of the recurrence law defining R% m , 

= E 2 m,o( a i4>) = dE2 " n ~ 1 ' 1 ( a , 4>) - ^-E2m-i,i{a,4>), 
dx\ ax\ 
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we use g^j- ^ 0, and remark^l S 2 m-i,-i(a, 0) = to obtain, 

#2ro-l,l( a > <W = 0- 
If we assume E2 m -i,2k+i(4'i °) — 0, substituting it in the recurrence law 

E 2m ,2k+i(a, 4>) = dE2 j Xi 1 ' 2k i a ' 4>) - 2 ( fc + l)^-£ ; 2m-i,2(fe+i)(a, 0), 
we obtain that 

^2m-l,2(fe+l)(Q!) 4>) = 0, 

and we conclude by finite induction. □ 

Corollary 3 Let H be a classical hamiltonian of type (|12[) . then the following 
statements are equivalent. 

1. The NVE for generic integral curve ()14fl lying on T has polynomial coef- 
ficient a(t) of degree 2m — 1. 

2. H is written, 

H = Vl — V2 + A - P2m-i(^i)^2 + P{xi,x 2 )xl 1 (35) 
for Ao constant, and f^m-lO^l) polynomial of degree 2m — 1. 

Proof. It is clear that condition 1. is satisfied if and only if (a, (j>) is a 
solution of i?2m and it is not a solution of i?2m-i- By the previous lemma, it 
implies J^L = o, and then the system Ri m is reduced to d dx 'i" and then, <f> is a 
constant and a is a polynomial of degree at most 2m — 1 . □ 

3.2.3 NVE Mathieu extended 

This is the standard Mathieu equation, 

£ = (co + ci cos(u;f))f , uj^O. (36) 

We can not apply our method to compute the family of hamiltonians corre- 
sponding to this equation, because {cq + c\ cos(u;i)} is not the general solution 
of any differential polynomial with constant coefficients. But, let us consider 

^, n d?a 9 da ,„„. 
Q(a) = w+ ^-, (37) 

the general solution of {Q(a) = 0} is 

a(t) = co + ci cos(wi) + C2 siri(ujt). 

Just notice that, 

c 2 



ci cos(uit) + C2 sin(oji) = \J c\+ c^ cos ^u>t + 



arctan 

ci 
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thus NVE l|14f) . when a is a solution of J37J), is reducible to Mathieu equation 
<|36l) by a translation of time. 

Using Proposition El we find the system of differential equations that deter- 
mine the family of hamiltonians, 

d 3 a da d 2 (j) ^ 2 da 

dx\ ' dx\ dx\ dx\ dx\ dx\ 

General solution of the first equation is 

a = A + \\X\ + X 2 x\. 

substituting it in the second equation, and writing y = we obtain a non 
homogeneous linear differential equation for y, 



dy 6X 2 y 



,2 



dx\ Ai + 2X 2 x\ 



= U! 



(38) 



We must distinguish two cases depending on the parameter. If X 2 = 0, then we 
just integrate the equation by trivial quadratures, obtaining 



2 2 

UJ xf 

Ho + /j,iX! H — 



and then, 



tt 2/l + 2/2 , i ^ Xi \ 2 \ 2 , of \ 3 ron\ 

H = h/i + Mi^i H 2 0X2 ~ *i x i x 2 + P(xi,x 2 )x 2 , (39) 

If 7^ 0, then we can reduce the equation to separable using 

6X 2 y 



obtaining 



Ai + 2X 2 xi ' 



3du 6X 2 dx 3uj 3/i 



and then 



iuo 2 -Au Ai+2A 2 a;i' 4 4(Ai +2A 2 xi) 4 ' 

1 / 2 2 Mi 

2/ = — w Ai + 2w A 2 xi + — — — 

8A 2 V (Ai + 2X 2 xi) A , 

and finally we integrate it to obtain 0, 

jjii I w 2 AiXi uj 2 x\ 



= / ydxi = no 



32A| (Ai + 2A 2 a;i) 2 8A 2 



scaling the parameters adequately we write down the general formula for the 
hamiltonian, 

2/1+2/2 , Mi , Aiw 2 a;i 

" = 7, r MO 



(A!+2A 2 2;i) 2 8A 2 8 
Aqx 2 — X\x\x\ — X 2 x\x\ + (3(xi, x 2 )x2- (40) 
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3.3 Application of non-integrability criteria 

3.3.1 NVE with a(t) polynomial 

According to theorem (see Appendix), the Galois group corresponding to 
equations, 

e = p(t)t 

where P(t), is a non-constant polynomial, is a connected non-abelian group. 
So that, we can apply theorem^] and we get the following result. 

Proposition 4 Hamiltonians, 

H = y±±yl + x 2 Q(x 1 ) + f3(x 1 ,x 2 )xj 

where Q(xi) is a non- constant polynomial, and /3(x%, x%) analytic function around 
V , do no admit any additional rational first integral. 

Corollary 4 Every integrable (by rational functions) polynomial potential with 
invariant plane F = {x 2 =2/2=0} is written in the following form 

V = Qi(xi, xz)x 2 + Xix 2 + Ao, Ao,AieC. 

3.3.2 NVE reducible to quantum harmonic oscillator 

We have seen, that hamiltonians (|33|l . has generic NVE along curves in T of 
type Once again, we apply theorem 

Proposition 5 Hamiltonians of the family do not admit any additional 

rational first integral. 

We can also discuss, the Picard-Vessiot integrability of those equations (|32|l . 
First, we shall notice that by just an scaling of t, 

t= T Cl 



we reduce it to an quantum harmonic oscillator equation, 



dr AJc% 



(r 2 - E)£, E = (41) 

4 

In appendix (remark [SJ we analyze this equation. It is Picard- Vessiot in- 
tegrable if and only if E is an odd positive number. Then, let us compute de 
parameter E associated to NVE of integral curves of Hamiltonians 1331 

Let us keep in mind that the family of those curves is parameterized by 

h= yf , a 4 



2 (X 2 +2X 3Xl y 
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In order to fix the parameterization of those curves, let us assume that time 
t = corresponds to x\ — 0. The NVE corresponding to a curve, depending on 
energy h, is written: 

|= (c {h) + Cl (h)t + c 2 (h)t 2 )£. 

We compute these coefficients Ci(h) using, 

_ v/2fe(A 2 + 2A 3 a;i) 2 -2A4 t-o \J1hX\ - 2A 4 
yi A 2 + 2A 3 a;i A 2 

and then, by applying the hamiltonian field, 

c o{h) = y, ci(fe) = \f~ ^ 2 2 Xi , c 2 (h)=X 3 h, 
and then, it is reducible to equation l|41|) with parameter, 
E = -JL_ ( X 2-4\ih _ A 4 \ 

8\/Al V v7? vW' 

If A 2 = 4A1A3 and A4 = 0, then parameter E vanish for every integral curve 
in r. For any other case, E is a non-constant analytical function of h. 

We have proven that those NVE are, generically not Picard- Vessiot integrable 
for any hamiltonian of the family. 

3.3.3 NVE Mathieu 

In order to apply theorem we just need to make some remarks on the field 
of coefficients. Let 7 be a generic integral curve of l|39[) . or Q4U|). Those curves 
are, in general, Riemann spheres. The field of coefficients -M 7 , is generated by 
x\,yx, so that it is C(xi,±i). We also have, 

a = A + Ai^i + X 2 Xi, a = i 1 (A 1 + 2A2X1), 

so that, for A2 = 0, we have 

M-y = (a, a) = C(smt,cost) = C(e lt ). 

and, for A 2 7^ 0, 

C{e lt ) <-+ M-y 

is an algebraic extension. 

So that, in our algebrization algorithm, the field of coefficients of Mathieu 
equation is taken, C(e lt ). For A 2 = 0, we can apply directly theorem^] and for 
A2 ^ 0, we can apply theorem El and then theorem 

Non trivial equations of type Mathieu, with field of coefficients C(e 1 *), ana- 
lyzed in remark [7| have Galois group SL(2, C). Thus for computed families of 
hamiltonians with NVE of type Mathieu, we get: 

Proposition 6 Hamiltonians of the families (|39|l and (|40|) if Ai ^ and 

(Ai,A 2 ) 7^ (0, 0) respectively, do not admit any additional rational first inte- 
gral. 
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Appendix 

A Kovacic's Algorithm 

This algorithm is devoted to solve the reduced linear differential equation (RLDE) 
£" = r£ and is based on the algebraic subgroups of SL(2, C). For more details 
see |Ko| . Improvements for this algorithm is given in |Ul-We| . in where is not 
necessary to reduce the equation. Here, we follows the original version given by 
Kovacic. 

Theorem 7 Let G be an algebraic subgroup of SX(2,C). Then one of the 
following four cases can occur. 

1. G is triangular izable. 

2. G is conjugate to a subgroup of infinite dihedral group ( also called meta- 
abelian group) and case 1 does not hold. 

3. Up to conjugation G is either of following finite groups: Tetrahedral group, 
Octahedral group or Icosahedral group, and cases 1 and 2 do not hold. 



Each case in the Kovacic algorithm is related with each one of the algebraic 
subgroups of SL(2, C) and the associated Riccatti equation 



According to theorem[7] there are four cases in the Kovacic algorithm. Only 
for cases 1, 2 and 3 we can solve the differential equation RLDE, but for the 
case 4 we have not Liouvillian solutions for RLDE. Is possible that the Kovacic 
algorithm only can provide us only one solution (yi), so that we can obtain the 
second solution (1/2) through 



I G = SL(2,C). 









(42) 
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Notations. For the equation RLDE with 

r = -, s,t e C[x] 

we use the following notations. 

1. Denote by V be the set of (finite) poles of r, V = {c e C : t(c) = 0}. 

2. Denote by T = V U {oo}. 

3. By the order of r at c e V , o(r) c , we mean the multiplicity of c. 

4. By the order of r at oo, o (roc) , we mean the order of oo as a zero of r. 
That is o (r^) = deg(t) — deg(s). 

5. By the order of r at c e V , o(r) c , we mean the multiplicity of c. 
A.l The four cases 

Case 1. In this case \y/r\ c and [\A*]oo means the Laurent series of y/r at 
c and the Laurent series of yfr at oo respectively. Furthermore, we define 
e(p) as follows: if p € T, then e (p) e {+,—}. Finally, the complex numbers 
a+ , , a+ , will be defined in the first step. If the differential equation has 
not poles it only can fall in this case. 

Step 1. Search for each c <G r" and for oo the corresponding situation such 
as follows: 

(co) If o (r c ) = 0, then 

[vf] c = 0, a±=0. 

(ci) If o (r c ) = 1, then 

[V^] c = 0, a± = l. 

(c 2 ) If o (r c ) = 2, and 

r = h b(x — c)~ 2 + • • • , then 

[\/r] c = 0, a± = . 

(c 3 ) If o (r c ) = 2u > 4, and 

r = (a (jc - c)~ v + ... + d(x- cy 2 ) 2 + b(x - c)-(" +1 ) + . . . , then 

[V^] c = a (»- c) - '" + ... + d{x- c)~ 2 , a± = ^ +w 
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(001) If o (roo) > 2, then 

(00 2 ) If o (roo) = 2, and r = • • • + 6a; 2 + • • • , then 

r n ± 1±a/T+46 

Knoo^ ' = 2 

(003) If o (roo) = — 2v < 0, and 

r = (ax" + ... + df + bx v - 1 + ■■ ■ , then 

IV^oo = ^ + ••• +C? ' and a ~ = ^( ± ^ _ ' y )" 

Step 2. Find L> ^ defined by 

D=^neZ + :m = a^°> - J>" (c \V( e (p)) per J . 

If D = 0, then we should start with the case 2. Now, if #D > 0, then for each 
raeflwe search w e C(x) such that 

W = e (00) [ Vr] x + ]T ( £ ( c ) W c + < (C) (x- tra- 
cer' 

Step 3. For each m <G D, search for a monic polynomial P m of degree m 
with 

P^ + 2uP' m + (w' + w 2 -r)P m = 0. 

If success is achieved then £1 = P m ef w is a solution of the differential equa- 
tion RLDE. If not, then case 1 cannot hold. 

Case 2. Search for each c <G T' and for 00 the corresponding situation such 
as follows: 

Step 1. Search for each c e V and 00 the sets E c ^ and 7^ 0. For 
each c e T' and for 00 it is define E c C Z and C Z as follows: 

(ci) Ifo(r c ) = l, then^ c = {4} 

(c 2 ) If o (r c ) = 2, and r = ■ ■ ■ + b(x — c)~ 2 + • • • , then 

E c = 1 2 + kVTTIb : k = 0, ±2} . 
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(c 3 ) If o (r c ) = v > 2, then £ c = {v} 

(001) If o (r-oo) > 2, then = {0, 2, 4} 

(00 2 ) If o (7*00) = 2, and r = ■ ■ ■ + bx 2 + • • • , then 

£oo = 1 2 + W 1 + 4& : fc = 0, ±2| . 

(003) If o (roo) = v < 2, then ^ = {v} 



Step 2. Find D ^ defined by 
D = < m £ Z + : m = 



5 ^ e °° - ^2 e cj e £ p , p e r j . 



If Z? = 0, then we should start the case 3. Now, if #D > 0, then for each m e D 
we search a rational function 9 defined by 



2^ x-c 
cer' 



Step 3. For each m £ D, search a monic polynomial P m of degree m, such 
that 

PZ + WP'Z + (36»' + 36» 2 - 4r)P r '„ + (<?" + 300' + 6» 3 - 4r<9 - 2r') P m = 0. 

If P m there is not exists, then the case 2 cannot hold. If such a polynomial is 
found, set <j> = 9 + P'/P and let lo be a solution of 

uj 2 + 4>oj + i (0' + </> 2 - 2r) = 0. 

Then £1 = ef " is a solution of the differential equation RLDE. 

Case 3. Search for each cer' and for 00 the corresponding situation such 
as follows: 

Step 1. Search for each cer' and 00 the sets E c 7^ and E^ ^ 0. For 
each ceT' and for 00 it is define E c C Z and E^ C Z as follows: 

(ci) If o (r c ) = 1, then E c = {12} 

(c 2 ) If o (r c ) = 2, and r = h b(x - c)~ 2 H , then 

E c = 1 6 + fev/T+lft : k = 0, ±1, ±2, ±3, ±4, ±5, ±6} . 



29 



(oo) If o (roo) = v > 2, and r = ■ ■ ■ + bx 2 + • • • , then 



{12k 1 
6 + — VTT45: fc = 0,±l,±2,±3,±4,±5,±6L ne {4,6,12}. 



Step 2. Find D ^ defined by 
D 



jm e Z+ : m = ^ ^eoo - ^e c ^ ,Ve p e E p , p e r| . 

In this case we start with n = 4 to obtain the solution, afterwards n = 6 and 
finally n = 12. If D = 0, then the differential equation has not Liouvillian 
solution because falls in the case 4. Now, if #D > 0, then for each to G D with 
its respective n, it is search a rational function 



12 ^ a; - c 
cer' 



and a polynomial 5 defined as 

5= n^- c )- 



cer' 



Step 3. Search for each raeD, with its respective n, a monic polynomial 
P m = P of degree m, such that its coefficients can be determined recursively by 

P_i =0, P„ = -P, 



- -5P/ - ((n - i) 5' - Sfl) P - (n - i) (i + 1) 5 2 rP +1 , 

where i € {0, 1 ... , n— 1, n}. If P there is not exists, then the differential equation 
has not Liouvillian solution because falls in the case 4. Now, if P there exists, 
it is search ui such that 

" qi p 

(«-*)' 

then a solution of the differential equation RLDE is given by 
where w is solution of the previous polynomial of degree n. 
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A. 2 Some remarks on Kovacic's algorithm 

Along of this section we assume that RLDE falls only in one of the four cases. 

Remark 12 (Case 1) If RLDE fall in case 1, then its Galois group is given 
by any of the following groups: 

II e when the algorithm provides two rational solutions or only one rational 
solution and the second solution obtained by equation Q42|l has not loga- 
rithmic term. 



this group is connected and abelian. 
12 Gfc when the algorithm provides only one algebraic solution £ such that 



this group is disconnected and its identity component is abelian. 
13 C* when the algorithm provides two non algebraic solutions. 



this group is connected and abelian. 

14 C + when the algorithm provides one rational solution and the second solu- 
tion is not algebraic. 



this group is connected and abelian. 

15 C* ix C + when the algorithm only provides one solution £ such that £ and 
its square are not rational functions. 



This group is connected and non-abelian. 

16 SL(2, C) if the algorithm do not provides any solution. This group is con- 
nected and non-abelian. 

Remark 13 (Case 2) If RLDE fall in case 2, then the Kovacic's Algorithm 
can provides one or two solutions. This depends of r such as follows: 




£ fe S C(ar) and ^ $ C(x). 
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Ill if r is given by 



r = 



2(f)' + 2cj>-cj> 2 
4 



then there exists only one solution, 



112 if r is given by 



2(f)' + 2cj>-(j> 2 
4 



then there exists two solutions. 
113 The identity component of the Galois group for this case is abelian. 

Remark 14 (Case 3) If RLDE falls in case 3, then its Galois group is given 
by any of the following groups: 

III1 Tetrahedral group when uj is obtained with n = 4. This group of order 
24 is generated by 



III2 Octahedral group when u> is obtained with n = 6. This group of order 
48 is generated by 



III3 Icosahedral group when uj is obtained with n = 12. This group of order 
120 is generated by 



III4 The identity component of the Galois group for this case is abelian. 
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